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In this article, we study the thermalizability of a system consisting of two atoms in a circular, 
transversely harmonic waveguide in the multimode regime. While showing some signatures of the 
quantum-chaotic behavior, the system fails to reach a thermal equilibrium in a relaxation from 
an initial state, even when the interaction between the atoms is infinitely strong. We relate this 
phenomenon to the previously addressed unattainability of a complete quantum chaos in the Seba 
billiard [P. Seba, Phys. Rev. Lett., 64, 1855 (1990)], and we conjecture the absence of a complete 
thermalization to be a generic property of integrable quantum systems perturbed by a non-integrable 
but well localized perturbation. 
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Introduction.- The ultracold quantum gases have been 
long proven to be an ideal testbench for studying the 
fundamental properties of quantum systems. One of the 
most interesting and frequently addressed topics is the 
behavior of quantum systems in the vicinity of an inte- 
grable point. Experimental results already include the 
suppression of relaxation of the momentum distribution 
[H and the modified decay of coherence 0. Note that 
the fully developed quantum chaos with ultracold atoms 
has been investigated in experiments 

Traditionally, the underlying integrable system is rep- 
resented by the Lieb-Liniger gas [4| that is well suited 
to describe the dynamics of Van-der-Waals interacting 
bosons in a monomode atomic guide (see [B[ for a re- 
view). The nontrivial integrals of motion there are inti- 
mately related to both the one-dimensional character of 
the atomic motion and to the effectively zero value of the 
two-body interaction range. The causes for lifting inte- 
grability include the virtual excitation of the transverse 
modes during the collision 0] and the coupling between 
two parallel weakly interacting Lieb-Liniger gases [2|. 

One would expect that in the multimode regime, when 
the motion becomes substantially multidimensional, no 
remnants of the former integrals of motion will survive. 
However, in this article we show that for the case of only 
two interacting atoms a strong separation between the 
transverse and longitudinal degrees of freedom remains, 
even in the case of the infinitely strong interaction be- 
tween the atoms. We attribute this effect to the short- 
range nature of the interaction. 

For our system we study both the degree of the eigen- 
state thermalization and the actual thermalization in an 
expansion from a class of realistic initial states. The 
eigenstate thermalization - the suppression of the 
eigenstate-by-eigenstate variance of quantum expecta- 
tion values of simple observables - provides an ultimate 



upper bound for a possible deviation of the relaxed value 
of an observable from its thermodynamical expectation, 
for any initial state in principle. However, recently a new 
direction of research has emerged: quantum quench in 
many-body interacting systems 0, 0, 0, E3] • In this class 
of problems, the initial state is inevitably decomposed 
into a large superposition of the eigenstates of the Hamil- 
tonian governing the dynamics, and the discrepancy be- 
tween the result of the relaxation and thermal values is 
expected to be diminished. In our article we show that 
in the case of two short-range-interacting atoms in a har- 
monic waveguide, a complete thermalization can never 
be reached under either scenario. 

We relate the absence of thermalization in our system 
to the previously addressed unattainability of a complete 
quantum chaos in the Seba billiard - flat two-dimensional 
rectangular billiard with a zero-range scatterer in the 
middle fllj j . Similarly to our system, the Seba billiard 
does show some signatures of the quantum-chaotic be- 
havior, although both the level statistics 12l4l4| and the 
momentum distributions in individual eigenstates 
show substantial deviations from the quantum chaos pre- 
dictions. Both the Seba billiard (first suggested and 
solved in Ref. [lH) and our system allow for an exact 
analytic solution, in spite of the absence of a complete 
set of integrals of motion. 

The system- Consider two short-range-interacting 
atoms in a circular, transversally harmonic waveguide. 
In this case the center-of-mass and the relative motion 
can be separated. Note that the harmonic approximation 
for the transverse confinement has been proven to model 
real life atom waveguides with high accuracy. Predictions 
of the analogous model involving an infinite waveguide 
15 , 17] has been confirmed by a direct experiment-theory 



comparison of the energies of the quasi-one-dimensional 
dimers 18| and numerical calculations involving anhar- 
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monic potentials (see review [f| for references). 

The unperturbed Hamiltonian of the relative motion 
is given by the sum of the longitudinal and transverse 
kinetic energies and the transverse trapping energy U — 
ticj 2 ± p 2 /2, 
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2fi dz 2 
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U — fiuj\ 
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where lu± is the transverse frequency, A p is the transverse 
two-dimensional Laplacian, and p is the reduced mass. 
In our model, the transverse and longitudinal degrees of 
freedom are coupled by a potential of a Fermi-Huang type 

1: 



V = 



2nfi 2 a s 



6 3 (r)(d/dr)(r-), 



(2) 



were, a s is the three-dimensional s-wave scattering 
length. This approximation is valid at low collision en- 
ergies (e.g., below 300 mK for Li, see a recent review 
19(, while the ultracold atom energies lie well below 
1 uK). The model ^ stands in a excellent agreement 
with the experimentally measured equation of state of 
the waveguide-trapped Bose gas [2(J and with numerical 
calculations with finite-range potentials 

The ring geometry of the waveguide imposes period- 
ic boundary conditions along z-direction. In what fol- 
lows, we will restrict the Hilbert space to the states that 
have zero z-component of the angular momentum and 
that are even under the z O — z reflection; the inter- 
action has no effect on the rest of the Hilbert space. 
The unperturbed spectrum is given by (see [HI]) E n \ = 
2Tuv±n + h 2 (2Trl/L) 2 /(2fi), where n > and I > are 
the transverse and longitudial quantum numbers, respec- 
tively. 

The interacting eigenfunctions with the eigenenergy 
E a can be expressed as (see [22I ]) 



(p,z\a) = C a 



:os {2\Je a — An£) 
n ^ ^e a -Xn sin ^/e a - Xn 



e-§4°)(0(,3) 



where the rescaled energy e Q is given by e a = 
XE a / (2huj±), X = (L/a±) 2 is the aspect ratio, £ = 

z/L— 1/2, £ = (p/a±) 2 , Ln (£) are the Legendre poly- 
1/2 

normals, a± — (fi/puj±) is the size of the transverse 
ground state, and C a is the normalization constant. The 
eigenenergies are solutions of the following transcenden- 
tal equation (see Q): 
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where £ [y,x) is the Hurwitz Q- function (see [Hj]). Similar 
solutions were obtained for two atoms with a zero-range 
interaction in a cylindrically-symmetric harmonic poten- 
tial [13] (that system was analysed numerically in Ref. 



24]). Higher partial wave scatterers were analyzed in 
the Ref. 

At rational values of X/tt 2 , the unperturbed energy 
spectrum shows degeneracies that are not fully lifted in 
the full [deduced from (jl])] spectrum. To minimize the 
effect of the degeneracies, we, following Ref. [Il|], fi x the 
length of the guide to (L/a±) 2 = X = 2(f> gr n 7 sa 9774, 
where 4> gr = (1 + Vo)/2 is the golden ratio. 

Standard quantum- chaotic tests- The results of the 
study of the level spacing distribution in our system are 
fully consistent with the analogous results for the Seba 
billiard 12, [l3| • For the energy range E > lOOTiaa 
and the aspect ratio (L/a±) 2 = 2(f) gr ir 7 , the distribu- 
tion quickly converges to the Seba distribution [13] at 
a s > 10 _1 a^. This distribution does show a gap at 
small level spacings but fails to reproduce the Gaussian 
tail predicted by the Gaussian Orthogonal Ensemble. At 
small a s the level spacing distribution tends to the Pois- 
son one. We have also verified that in the unitary regime, 
a s 3> ax, the statistics of the point-by-point variations 
of the eigenstate wavefunction in the spatial representa- 
tion is close to the Gaussian one, according to the general 
prediction [2H ] and a particular observation in the case 
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of the Seba billiard 

Eigenstate thermalization— According to the eigen- 
state thermalization hypothesis 0-H, the ability of an 
isolated quantum system (with few or many degrees of 
freedom regardless) to thermalize follows from the sup- 
pression of the eigenstate-by-eigenstate fluctuations of 
the quantum expectation values of relevant observables. 
Scattered dots at the Fig. [1] show the quantum expec- 
tation values of the transverse trapping energy, (a|[/|a) 
as a function of the eigenstate energy. The variation 
of (q|£/|a) remains (in contrast to the quantum-chaotic 
billiards [27[) of the order of the mean, even for the en- 
ergies much larger than any conceivable energy scale of 
the system. Here and below, we work in the regime of 
the infinitely strong interactions, a s — 10 6 aj_. 

For a given eigenstate, quantum expectation values 
of the occupation probabilities for the transverse (Fig. 
Hk) and longitudinal (Fig. [2Jd) modes show no conver- 
gence to their microcanonical expectation values, P(n) oc 
/ dl 5(E n i — E) and P(l) cx Jdn 5(E„i ~ E) respectively. 
Both distributions are heavily dominated by only a few 
peaks (cf. scars in the momentum space for Seba billiards 

El)- 

Relaxation from an initial state - Now we are going 
to address directly the ability of our system to ther- 
malize from an initial state \ip(t = 0)). In this case, 
the infinite time average of the quantum expectation 
value of an observable A will be given by A rc i ax = 

hm^^a/w,) r o 4 — (ip(t)\A\m) = £J<c#(* = 

0)| 2 (o|A|n:), where \a) are the eigenstates of the system 
(see [8(|). For a quantum-chaotic system, it is expected 
that the above infinite time average coincides with the 
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FIG. 1. The time average of the transverse trapping energy 

^rclax. = ^o.t maK -Kx>0-/tmax) J max (ip(t) a f ter relax- 
ation from an initial state as a function of the energy of the 
initial state E [see Hamiltonian ([I}]. Four families of the ini- 
tial states are considered. Long-dashed (green) line: state © 
with no = 0, 8 = 0.99, and the energy being controlled via 
scanning lo; short-dashed (blue) line: the same as the previ- 
ous one, except for 8 — .1; dotted (purple) line: the state ([5]) 
with 8 = 0.1 and lo — and the energy control via no; solid 
(red) line: state ((6J with K2 = 2ki, 8 — 0.99, lo = and the 
energy control via ki. Legends schematically show the distri- 
bution of the initial states over the quantum numbers; arrows 
indicate the direction in which the parameters controlling the 
energy are scanned. The asterisk (*) corresponds to the set 
of parameters used to produce the Fig. [2] Grey dots show 
the quantum expectation value of the trapping energy U for 
every hundreds eigenstate of the full Hamiltonian (1 1121) . All 
the values of the trapping energy are shown with respect to 
its microcanonical expectation value Utherm. ~ E/3. For all 
the points, (L/a±) 2 = 2<j) gr n 7 « 9774 and a s = 10 6 a±. 



thermal prediction. 

Consider the following initial state: 



(z\ij)(t = G)) oc cos 



|C| cos(27rfoC)|no). (5) 



Longitudinally, the state is represented by the ground 
state of a length L5 hard-wall box split initially by an 
ideal beamsplitter with momenta ±2t:Iq/L; the box is 
centered at the maximal interatomic distance. The state 
is distributed among approximately tv/S axial modes lo- 
calized about I = ±Iq. The initial transverse state is 
limited to a single mode no. Note that for no = and 
5 « 1 the state is conceptually similar to the initial state 
used in the equilibration experiments 

After relaxation, both transverse (Fig. (3l) and longi- 
tudinal (Fig. [2}}) distributions remain very far from the 
equilibrium predictions. Fig. [TJ shows the equilibrium 
value of the transverse trapping energy U (see ((TJ) for 
three sequences of the initial states of the type ((5j. The 
first two, with low transverse energy, show an approx- 
imately -40% deviation of the relaxed value of U from 
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FIG. 2. Distribution of the mean occupation numbers for 
the transverse (a) and longitudinal (b) modes of the unper- 
turbed Hamiltonian. The long-dashed (green) and solid (red) 
lines show the initial and time-averaged distributions for an 
initial state © with no = 0, 8 = 0.99, and l = 292.8. En- 
ergy of the initial state is E = (2/ X)hu±e a = 173.2huj±. The 
short-dashed (blue) line corresponds to the 17096's eigenstate 
(closest in energy to E from the above, with a relative dis- 
crepancy of 10 _J ) of the interacting Hamiltonian. The dotted 
(purple) line shows the microcanonical prediction. The rest 
of the parameters is the same as at the Fig. [TJ 



the thermal prediction, thus being seemingly correlated 
with its initial, negative value. Note that the correla- 
tion between the final and initial values of observables 
in quantum systems is addressed in Ref. [2^]. Also, note 
that the initial states of the second sequence have a much 
greater spread over the longitudinal modes than the first 
one; consistently, the energy-to-energy variation of the 
relaxed values of U is less than for the first sequence. The 
third sequence, with low longitudial energy, shows devi- 
ations from the thermal prediction that range between 
+50% and +5% but never reach zero. 
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Another type of the initial state 

(p, 2#(t=0)) « cos ^9 f | - |C| ) r>» r>~U ) 



x e 



-Kif _ -re 2 f 
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(6) 



We are grateful to Felix Werner for the enlightening 
discussions on the subject. This work was supported by 
a grant from the Office of Naval Research (N00014-09-1- 
0502). 



allows for an additional spread over the transverse modes. 
The transverse wavefunction vanishes at the waveguide 
axis. The mode occupation has a minimum at n = 0, and 
the occupation of the ground transverse mode tends to 
zero for Ki < K% -C 1. The results for a sequence similar 
to the third sequence described above are shown at Fig. 
[T] as well. In spite of the significant initial spread over 
both transverse and longitudinal modes, the deviations 
from the equilibrium are close to the ones for the third 
sequence, while the energy-to-energy variations for the 
former are indeed less than for the latter. 

Note that we have also successfully tested the conver- 
gence of the instantaneous quantum expectation values 
of the observables to their infinite-time averages. 

Summary and interpretation of results - In this arti- 
cle, we solve analytically the problem of two short-range- 
interacting atoms in a circular, transversely harmonic 
multimode waveguide. We assess the ability of the sys- 
tem to thermalize from an initially excited state; a broad 
class of initial states has been analyzed. We find a sub- 
stantial suppression of thermalization, even for the in- 
finitely strong interactions. We associate this effect with 
the previously demonstrated unattainability of complete 
quantum chaos in the Seba- type billiards [llMl5j|. 

We conjecture that the effect of suppression of thermal- 
ization is generic for the integrable systems perturbed, 
no matter how strong, by a well localized perturbation. 
The following reasoning applies. In a quantum chaotic 
system, a given eigenstate \a) consists of a large superpo- 
sition of the eigenstates \n) of the underlying integrable 
system, which are drawn indiscriminately from the micro- 
canonical shell. Number of principal components in such 
a superposition is a sensitive measure of the approach to a 
complete chaos and the subsequent thermalizability (see, 
for example [29|). The zero-point-localized perturbation, 
being a particular case of a separable perturbation, gen- 
erates the eigenstates of a form (n\a) oc l/(E a — Eft) 
(see [22j). There is always one and only one perturbed 
ener gy E a in between any two unperturbed energies En 
(see [13|]). In addition, for the case of the infinitely strong 
perturbation (where the approach to chaos is expected to 
be the closest) E a tends to the middle position between 
the two Eft's. In this case, the only energy scale that 
remains is the distance between the unperturbed levels, 
and, as a result, the number of principal components a 
given eigenstate \a) consists of is always of the order of 
unity. This property extends to the case of finite-range 
interactions, whenever the interaction range is much less 
then the de Broglie wavelength of the colliding atoms. 

Our study is the first attempt to address thermaliz- 
ability of a quantum system with separable interactions. 
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The waveguide system as well as the Seba billiard 
0|, belongs to the class of problems where an integrable 
system of Hamiltonian Hq is perturbed by a separable 
rank I interaction V\£)(1Z\ (see e.g. Q). Eigenstates \a) 
of the interacting system are solutions of the Schrodinger 
equation 



[H + V\£)(R\}\a) =E a \a) 



(1) 



where E a are the corresponding eigenenergies. Let us 
expand the eigenstate \a) over the eigenstates \n) (of en- 
ergy Eft) of the non- interacting hamiltonian Hq. The 
Schrodinger equation leads to the following equations for 
the expansion coefficients 



(E a - E n )(n\a) = (n\£)(TZ\a) cx (n\£), 



(2) 



where the omitted factor in the second equality is inde- 
pendent of n. 

As a result, any eigenfunction \a) functionally co- 
incides with the Green function of the non-interacting 
Hamiltonian taken at the energy of the eigenstate E a , 



E 



\n)(n\£) 



E n 



Efi 



(3) 



The omitted factor is determined by the normalization 
conditions. 

Substitution of this expression to the Schrodinger 
equation dXJ) leads to the following eigenenergy equation: 



E 



(TZ\n)(n\£) 



E n 



Ea 



1 

V' 



(4) 



Similar expressions were obtained in Refs. 0, 0| for the 
case of the Seba billiard and its generalizations. 

In the case of the relative motion of two short-range- 
interacting atoms in a circular, transversely harmonic 
waveguide, the derivation uses ideas of the analogous 
model involving an infinite waveguide 0-0] • The unper- 
turbed Hamiltonian here is given by Eq. (1) in [J and the 



Fermi-Huang interaction (see Eq. (2) in can be writ- 
ten in the separable form [3j with the interaction strength 
and formfactors given by 



V = 



2nh 2 a<. 



\£)=S 3 (r), (ll\=5 3 (r)~(r-). (5) 



In what follows, we will restrict the Hilbert space to 
the states of zero z-component of the angular mo- 
mentum and even under the z o — z reflection; the 
interaction has no effect on the rest of the Hilbert 
space. The non-interacting eigenstates are products 
of the transverse and longitudial wavef unctions. The 
transverse two-dimensional zero-angular-momentum har- 
monic wavef unctions, 



1 



(p|n) = -^£Cp) (0exp(-C/2). 



(6) 



labeled by the quantum number n > 0, are expressed 
in terms of the Legendre polynomials, Ln ^ (£), where 
£ = (p/a±) 2 and ai = (h/ fiu^) 1 ^ 2 is the transverse os- 
cillator range. The longitudial wavefunctions, labeled by 
I > 0, are the symmetric plane waves satisfying periodic 
conditions with the period L: 

(z\k) = (2/L) 1/2 cos 27rfcC, (*|0) = L- 1 ' 2 , (7) 



where £ = z/L — 1/2. 
therefore given by 



The unperturbed spectrum is 



E nt = 2hco ± n + h 2 {2irl/L) 2 /(2fi) 



(8) 



Substituting the above non-interacting eigenstates and 
eigenenergies to Eq. (0} and using Eq. (1.445.8) in [§] 
for summation over I, one obtains Eq. (3) in [l[ for the 
interacting eigenstates. The normalization factor C a is 
determined by the condition 



2?r 



p 2 dp(a'\p, z)(p, z\a) = S ac 



(9) 
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Using orthogonality of the Legendre polynomials it can 
be expressed as 



1 



/ COt via"" ^ n 
-1/2 



E 



(e Q — An) sin y/e a — Ar 



(10) 



where the rescaled energy e Q is given by e a = 
A-Bq,/ (2?iu;_i_) and A = (L/a±) 2 is the aspect ratio. 

For the Fermi-Huang interaction the eigenenergy equa- 
tion ((H) attains the form 



d_ 

Or 



•E 



(p|n)(z|Q(n|0)(0|0) 

E — Eni 



1 

V 



(11) 



r=0 



In the limit r — > 0, the eigenstate is spherically symmet- 
ric. This allows us to deal with the z-axis only (see 0]). 
Substitution of the nonintcracting wavefunctions ((6|) and 
([7]) and energies (jSJ) with the subsequent summation over 
I leads to 



Va 



d_ 

dz 



n=0 



■S (2Ve a - An(*/L - 1/2)) 
y/e a — An sin y/e a — Xn 



z=0 



a± ' 
(12) 



This sum contains both the regular part and the irregular 
one, the latter being proportional to z _1 . The regular 
part can be extracted using the identity 



lim 

z->0 



^(An - e)- 1 ' 2 exp(-2 v / An - ez/L) 



,n=0 



(13) 



(see [6fl) where C [y, a) is the Hurwitz zeta function (see, 
e. g., [3,Q. Finally we arrive at the transcendental equa- 
tion for the eigenenergies (4) in [![. The summands in 
the sums Eqs. (3) and (4) in Q and in Eq. (fT0|) decay 
exponentially with n, leading to the fast converging se- 
ries. Note that the imaginary parts of the two terms 
in the left hand side of Eq. (4) in [l[ cancel each other 
automatically. 
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